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Abstract: The molecular motions of Fl-ATPases derived from bovine heart mitochondria (MF1) and thermoalkaliphilic 
bacteria (TA2F1) were studied in this paper using Gaussian Network Model (GNM)[13-15]. The data obtained fromGNM 
calculation using a cut off distance 10 A were compared and analyzed. MF1 was taken as the reference protein for 
comparison. The calculated force constants (y) for TA2F1 and MF1 proteins were 0.054 kcal/ (mole ■ A 2 ) and 0.08 
kcal/(mole- A 2 ) respectively. The larger force constant of MF1 implies that the inter residue contact potential of MF1 is 
stronger than that of TA2F1. The fluctuations exhibited by the different chains of aand /? subunits of MF1 are very high 
compared to the TA2F1 protein. The largest fluctuations ~ 2.5A 2 was observed in the carboxy terminal region of the chain B 
of MF1 suggesting that the a 3 B 3 cylinder of MF1 is flexible whereas the a 3 /? 3 cylinder of TA2F1 is highly rigid. The 
ysubunit of TA2F1 generated a fluctuation ~ 2. 5 A 2 . The fluctuation of esubunit of TA2F1 tagged with the protruding portion 
of the ysubunit is the highest ~3.5A 2 compared to all other subunits of both the proteins. The theoretical and experimental B 
factors match well. The correlation coefficient for MF1 was 0.71 (71%) and for TA2F1 was 0.64 (64%). In the case of both 
proteins, the hydrophobic 'sleeve' forming the molecular bearing was found to be located in the minima region of fluctuation 
curves. 
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I. INTRODUCTION 

Molecular motor proteins are nanoscale natural machines found in biological cells. They convert chemical energy 
derived from the hydrolysis of adenosine triphosphate (ATP) into mechanical workor motion that forms the basis of cellular 
function. ATP molecules serve as remarkable sources of chemical energy that is used to perform the various cellular 
functions. Scientists of diverse background are studying the structure and function of these motors to either to mimic their 
action or to incorporate them into novel devices. Motor proteins are involved in ATP synthesis, organelle transport, muscle 
contraction, protein folding and nucleic acid translocation. They also take part in cell signaling, division and motion. The 
most prominent examples of ATP-driven linear motors are myosin, kinesin and dynein. Myosin is involved mainly in muscle 
contraction and it slides on actin filament. Kinesin moves along microtubules and is responsible for transporting cellular 
cargo such as organelles and signaling molecules. Dynein is another microtubule based motor protein that does many cellular 
functions. Examples of rotary molecular motors are bacterial flagella motor and ATP synthase. 

The rotary enzyme ATP synthase found in mitochondria, bacteria and chloroplasts generates ATP using the free energy 
stored in transmembrane ion gradients [1-6]. This enzyme consists of two parts (F 0 and F[). The asymmetric membrane- 
bound F 0 that contains the proton channel and the soluble F[ that contain three catalytic sites which work cooperatively. 
These two motors are connected through an asymmetric rotor shaft made of the coiled-coil gamma subunit. The proton flow 
through F 0 leads to the ATP synthesis in Fj and ATP hydrolysis in Fj leads to the rotation of the rotor in opposite direction 
that finally pumps the protons back towards F 0 . The Fj portion consists of a 3 p 3 ye5 subunits. F, only hydrolyzes ATP and is 
called Fj-ATPase. Both the three a and three Psubunits can bind ATP molecules but only three P subunits contain catalytic 
sites in which ATP can be hydrolyzed [7]. The e subunit acts as an inhibitor of ATP hydrolysis in F r ATPases [8]. The 5 
subunit is a part of stator in ATP synthase [9]. ATP hydrolysis driven rotation of y subunit around 013P3 cylinder was 
observed in thermophilic FpATPase [10]. ATP was generated in P subunits of F r ATPase by mechanical rotation of the y 
subunit in a direction opposite to the direction of ATP hydrolysis [11]. 

The proteins are flexible and can adopt different conformations in response to their environmental conditions. The 
protein function is a dynamic property [12]. To elucidate protein function it is necessary to understanding structure and 
dynamics of proteins. The structures of different proteins in the folded or native state seem to be in static form but proteins 
are engaged in functional motions during various cellular activities. Such motion in proteins concertedly occurs within amino 
acids, subunits and domains. A protein in motion frequently changes its conformation to do a particular function. Many 
biological processes like catalysis, transport of metabolites, folding, binding, molecular recognition etc. are performed by 
different proteins in the presence of motion. 

Highly efficient as well as physically realistic computational methods and tools are essential to understand the 
fluctuation dynamics of proteins. The systematic computational characterization of protein dynamics by molecular dynamics 
simulations is highly time consuming. Moreover, in the case of all atom models, computation is very expensive due to the 
use of complicated potential energy functions. In order to overcome such problems, a very simplified Gaussian Network 
Model (GNM) was developed [13-15] to predict the fluctuations of protein structures in coarse-grained level. This model 
involves normal mode analysis (NMA) to assess collective motions of proteins. The GNM was constructed based on inter- 
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residue contact topology in the folded state with a single parameter potential. The observed motion in a protein is dissected 
into a collection of modes which give information on the molecular mechanisms relevant to biological function [16, 17]. 

There are mainly two types of modes in the case of protein motion: slow mode and fast mode. In a slow (low 
frequency) mode, the different parts of a protein undergo motions with large-amplitudes. The slow mode motions are also 
called 'global motions' that are essential for understanding the protein function [18-19]. The slow modes make major 
contributions to thermal conformational fluctuations [20, 21]. Such motions can influence the interactions of proteins with 
other molecules and its environment. The fast modes with higher frequency fluctuations are more localized that include only 
a few residues and can play an important role in signal transmission, enzyme reactions and other internal processes [18]. In 
some studies [13, 22], it was revealed that the minima in the global mode shapes coincide with the hinge sites of the protein 
molecule and the maxima correspond to substrate recognition sites. 

Different experimental methods are used to determine the structures of proteins such as X-ray diffraction, nuclear 
magnetic resonance (NMR), electron microscopy etc. The X-ray crystallographic experiments give useful structural 
information about thermal and other fluctuations of the atoms in a protein. In X-ray crystallographic method, the average 
position of each atom is determined based on the electron density map and every atom is assigned to have a B-factor whose 
magnitude is directly proportional to the mean square displacement from its mean position. The B-factors reflect the 
flexibility and dynamics of protein structures. Models for protein dynamics provide insights into the mechanisms of motion 
of large proteins [12]. The B-factors measured by experiments can give limited information about molecular collective 
motions relevant to biological functions [23]. 

The aim of the present work is to study the fluctuation dynamics of thermoalkaliphilic F r ATPase (denoted as TA2Fj) 
[24] by the GNM and compare the results with the data obtained for the reference protein mitochondrial F r ATPase (denoted 
as MFi) [7]. The structure files for MF, (PDB ID: lbmf) and for TA2F; (PDB ID: 2QE7) proteins were downloaded from 
Protein Data Bank (PDB) [25]. The mean square equilibrium fluctuations as well as B factors were calculated using oGNM 
server maintained by Yang et al. [26, thanks to Prof. I. Bahar, University of Pittsburgh, USA for her kind permission to use 
the server]. The calculated B factors and fluctuations for both proteins were analyzed and compared. It was found that the 
calculated values of B factors agree well with the experimental ones. 

The plan of this paper is given as follows: Section 2 deals with the Comparative Analysis of Data and section 3 is 
presented to discuss the work. Section 4 is devoted to the conclusion of the work. 

II. DATA ANALYSIS 

The mean square fluctuations and B factors of thermoalkaliphilic F r ATPase (TA2F[) and mitochondrial F r ATPase 
(MF,) were calculated using oGNM server [26]. The data obtained by calculation was thoroughly analyzed. 
2.1 Analysis of mean square fluctuations 

The mean square fluctuations exhibited by the different subunits of both the MF[ and TA2F! proteins were compared 
and analyzed in the following sections. 

2.1.1 Analysis of the mean square fluctuations for different chains of the a subunit 

Proteins take part in various cellular processes to keep the cells alive. The functions of the proteins such as catalytic 
activity, folding, binding and molecular recognition involve protein motions or fluctuations [8]. Determination and 
understanding the three dimensional structures of proteins in the folded forms is necessary because their function depend on 
their folded native forms. The motions of proteins or conformational changes relevant to biological events are called 
functional motions or functional dynamics [12]. The fluctuations involved in biological functions of proteins are generated in 
collective manner that consists of correlated/anti-correlated motions of large groups of atoms, residues or domains of the 
proteins. Studying fluctuation dynamics of proteins is thus important to know their functional mechanisms. 

The Gaussian Network Model (GNM) [15, 26] is used in this thesis to theoretically calculate the mean square 
fluctuations and B factors of the residues of two proteins: thermoalkaliphilic F r ATPase (TA2F,) and mitochondrial F r 
ATPase (MF,). The GNM dissects protein motions into a collection of various normal modes and the slowest modes usually 
give information on the collective motions relevant to a biological function [12]. The data obtained for these two proteins 
from GNM calculation [26] for the l sl slowest mode of motion were analyzed. A cutoff distance r c = 10 A for inter residue 
interactions was used for both the proteins. It was found in a previous study that GNM results are insensitive to the cutoff 
distance provided in the range 7.3 A<r c <15 A [12]. By GNM calculation, a force constant of y= 0.054 kcal/(mole-A 2 ) was 
obtained for TA2Fj and y= 0.08 kcal/(mole-A 2 ) was for MF[. The force constant is a measure of the strength of 
intramolecular potentials that stabilize the native fold [7]. A force constant of 1.0 + 0.5 kcal/(mol-A 2 ) has been obtained in 
some proteins in a previous investigation [29]. 

The structure of thermoalkaliphilic F r ATPase (TA2F0 consists of eight chains (A, B, C, D, E, F, G and H) and 
mitochondrial F r ATPase has seven chains (A, B, C, D, E, F and G) as discussed inthestructure of F r ATPase. The motion of 
the amino acid residues contained in different subunits of the two proteins is revealed in the fluctuation curves. In order to 
compare, the mean square fluctuations of the amino acid residues involved in three a subunits (A, B and C chains), three [3 
subunits (D, E and F chains) and one y subunit (G chain) of TA2Fi and MFi proteins were plotted against the relevant 
residue indices as shown in Figs. 2.1-2.3. All these curves were plotted on the same scale of the mean square fluctuations and 
residue indices. The motion of the amino acid residues contained in the e subunit (H chain) of TA2Fi protein is shown in Fig. 
2.4. 

The fluctuations of the three a subunits (A, B and C chains) of both the TA2F[ and MFi proteins are shown in Fig. 2.1 
in which there are three curves for each chain. In the top figure obtained for the chain A, the red curve represents the 
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fluctuations for TA2F] protein. This curve is almost straight and parallel to the residue index axis. The black line was 
obtained for MFi protein and it contains some significant fluctuations at different residue index positions. Comparing these 
two curves (red and black) on the same scale, it is clear that the fluctuation for the TA2Fi protein is almost 0. 5A 2 less than 
that of the TA2Fi protein around the carboxy terminal region. The amino terminal of TA2Fj protein starts at 27 (GLU) and 
ends at 500 (PRO) in the carboxy terminal region and for the MFj protein, the chain starts at 24 (ASP) and ends at 510 
(ALA). The green curve is an expanded version of the red curve for TA2F[ protein. The purpose of expansion is to observe 
whether there is any fluctuation in the red curve. As seen, this curve (green) contains almost the same number of maxima and 
minima as the black curve. The positions of the maxima and minima also coincide with the black curve. That is, the pattern 
of both curves is almost the same. 

The middle and bottom figures contained in Fig. 2.1 were plotted for the chains B and C of a subunits. Similar to the 
top figure for the chain A, these two figures also possess three curves. In the middle figure for the chain B, the curve (red) 
for the fluctuations of TA2Fj protein is parallel to the residue index axis except around the 400 position in which a peak 
appears. In the expanded version (green curve), this peak appears to be the largest one among all other maxima. On the same 
scale of the red curve obtained for the TA2Fi protein, the black curve for the reference protein MFj contains many maxima 
and minima, the largest fluctuations occur after the 400 position of the residue index axis. The value of the fluctuations for 
MFj protein (black curve) is always higher from amino terminal to the carboxy terminal regions of the chain B (residues 
0-500) compared to the fluctuations obtained for the TA2F[ protein except near 300 position where a minima of zero 
fluctuations appears that coincide with the red curve. The largest variation of the fluctuations in the red and black curves of 
these two proteins is about 2. 5A 2 near the 500 position in the carboxy terminal region of the chain B. The amino terminal of 
the B chain for the TA2F t protein started at 27 and ended at 500 in the carboxy terminal whereas the amino terminal in the 
same chain for the reference MF, protein started at 24 and ended at 510 in the carboxy terminal region. Thus the number of 
amino acids in the B chain for the TA2F] protein and MF! protein differs by a factor of 13 and the amino terminal of the B 
chain of the reference MF, protein leads by 3 amino acid residues compared to the amino terminal of the B chain of the 
TA2F] protein. It seems to be that the maxima and the minima contained in the B chain for both the proteins do not coincide 
rather they reside in opposite positions, although in the case of chain A, the maxima and minima for both the proteins 
coincide. However, if the curve for the TA2Fi protein (say green) is moved toward the (amino terminal end) left on the 
residue axis, the trend of maxima and minima is almost similar to that of the reference protein MF! up to around the 400 
position. After this position, the trend is totally different in the carboxy terminal region which means that the black curve 
(MFi) although contains several minima but it shows its highest level of fluctuations. In the case of TA2F) protein, this 
region is almost flat exhibiting zero fluctuations except a little maximum at around the 435 position. 

The red curve (27GLU-500PRO) in the bottom figure obtained for the chain C of a subunit represents the fluctuations 
of TA2F] protein. The fluctuation is very little throughout the curve when compared with the black curve (19ALA-510ALA) 
containing a fluctuation range 0-1 .oA 2 " plotted for the MF! protein on the same scale. In the expanded curve (green), it is 
clear that some maxima and minima do not coincide but both the curves (black and green) show the similar trend. 
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Figure 2.1: The mean-square fluctuations exhibited by different chains of a subunit due to the 1 slowest mode were plotted 
as a function of residue index. The mean-square fluctuations were multiplied by 1000. The top, middle and bottom curves are 
for the chains A, B and C respectively. The red curve is for the TA2F] protein, the black curve is for the MF] protein and the 
green curve is an expanded version of the red curve. The fluctuations in right axis are for the green curve. The data were 
plotted using Sigma Plot. 

2.1.2 Analysis of the mean square fluctuations for different chains of the P subunit 

The chains (D, E and F) of the P subunit also undergo fluctuations as shown in Fig. 2. 2. The top figure contains three 
curves. It appears that the red curve (2ASN-462LEU) obtained by plotting the data for the D chain of TA2F] protein shows 
very little but constant fluctuations compared to the black curve (9THR-475GLU) for the MF] protein. This contains ~ 0.5 
A 2 fluctuations higher than the red curve. The green curve is an expansion of the red curve which shows similar fluctuations 
as the black curve. A maximum of green curve near 400 residue position seems not to be similar to the black curve because a 
minimum is observed in the black curve around this residue position. 

The middle figure is for the E chain. The red curve (2ASN-462LEU) for the TA2F] protein exhibits almost the same 
fluctuations up to 360 residue position compared to the black curve (9THR-474ALA) for the MF] protein. A peak appears at 
383 residue position and then the curve is flat up to the end. From the start to the end, the red curve contains fluctuations in 
the range of 0-0.75A 2 . The black curve have 0-0.5A 2 fluctuations and involves some maxima and minima throughout the 
length of the curve. In the expanded green curve, the maxima and minima contained in the red curve are shown. When the 
green and red curves are compared, it appears that the fluctuation pattern of the E chain of TA2F] protein differs from MF! 
protein because the positions of the maxima and minima are not the same. 

The three curves in the bottom figure were plotted for the F chain of the P subunit. The black curve (10THR- 
474ALA) for the MF] protein have fluctuations in the range of 0-0.46A 2 . The red curve (2ASN-462LEU) for the TA2F] 
protein is parallel to the residue axis. The green curve shows the fluctuations contained in the red curve. Comparing the 
black and green curves, it is evident that the fluctuation pattern in the F chain of both the proteins is similar. 




Figure 2.2: The mean-square fluctuations exhibited by different chains of P subunit due to the 1 st slowest mode were plotted 
as a function of residue index. The mean-square fluctuations were multiplied by 1000. The top, middle and bottom curves 
are for the chains D, E and F respectively. The red curve is for the TA2F] protein, the black curve is for the MF] protein and 
the green curve is an expanded version of the red curve. The fluctuations in right axis are for the green curve. The data were 
plotted using Sigma Plot. 

2.1.3 Analysis of the mean square fluctuations for the y subunit 

The y subunit of both the TA2Fi and MFi proteins contain only one chain which is denoted by the letter G. The 
fluctuations exhibited by the G chain are displayed in Fig. 2. 3. The red curve (3GLY-266ALA) obtained for the TA2F] 
protein shows fluctuations in the range of 0~2. 7A 2 which is much higher than the fluctuations exhibited by the six chains (A, 
B, C and D, E, F) of a and P subunits (Figs.2. 1-2.2). The black curve (1ALA-272LEU) is obtained for the G chain of MF, 
protein and the cyan curve is an expansion of the black curve. In the cyan curve, there is a very little fluctuation in the amino 
terminal around the 1-20 residue positions and at the end of the carboxy terminal of the G chain which is insignificant when 
compared with the fluctuations in the red curve. It is clear from the figure that the motion of the amino acid residues of G 
chain of TA2F] protein is far different from that of MF] protein. 
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Figure 2.3: The mean-square fluctuations as a function of residue index, driven by the 1 st slowest mode for y subunit. The 
mean-square fluctuations were multiplied by 1000. The red curve is for the TA2F, protein, the black curve is for the MF, 
protein and the cyan curve is an expanded version of the black curve. The fluctuations in right axis are for the green curve. 
The data were plotted using Sigma Plot. 

2.1.4 Analysis of the mean square fluctuations for the £ subunit 

The structure of TA2F] protein contains the subunit e having one chain which is expressed by the letter H.The 
fluctuation of this subunit is plotted in red color as shown in Fig. 2. 4 and it has a fluctuation of ~3. 7A 2 which is the highest 
level of fluctuation among all the chains (A, B, C, D, E, F and G) of other subunits (a 3 p 3 ). 




Figure 2.4: The mean-square fluctuations of the £ subunit (H chain) of TA2Fi protein were plotted as a function of residue 
index. The protein motion was driven by the 1 st slowest mode. The mean-square fluctuations were multiplied by 1000. The 
data were plotted using Sigma Plot. 



2.2 Analysis of B factors 

The GNM theory was used to calculate the B factors or Debye-Waller factors or temperature factors of individual 
residues of different chains of MFi and TA2Fj proteins. The calculated B factors were compared with the experimentally 
measured B factors as reported in the PDB files (PDB ID:lbmf for MFi and PDB ID:2QE7 for TA2F,) as shown in Figs.2.5- 
2.12. The B-factors are related to the mean square fluctuations of individual residues according to the equation 4.10 in 
chapter IV. The calculated correlation coefficient for MFi was 0.71 (71%) and for TA2F! was 0.64 (64%). The correlation 
coefficient for ribonuclease Tl was calculated as -0.7 (70%) by GNM theory in a previous work [30]. 

In general, the overall agreement of calculated B factors with the experimental values for the chains A, B, C, D, E 
and F of both the MFj and TA2F, proteins is very good although there are some deviations in calculated maxima and minima 
from experimental ones for individual residues (Figs. 2. 5-2. 10). Some examples are given here. In the case of chain A (blue 
and pink curves) of MF, protein, there are some deviation of theoretical maxima around 222 ASP, 315ALA, 368GLY and 
510ALA and minima around 323ALA and 403PHE from experimental ones. For the same chain (dark cyan and dark yellow 
curves) of TA2F, protein, significantly deviated maxima and minima are 57ASN, 69ASP, 310GLY and 456PHE, 482PRO, 
492GLU respectively. For other chains (B, C, D, E and F) of MFi protein, it is seen from the top figures that maxima: 
222ASP, 314ASP, 360GLY and 506ALA (chain B), 20 ASP, 58GLY, 222ASP, 314ASP and 510ALA (chain C), 41ARG, 
208LEU, 361ASN and 450ASP (chain D), 40GLY, 209LYS and 473LEU (chain E) and 40GLY, 209LYS and 360PRO 
(chain F) deviate significantly. Minima: 35GLY and 341 ASN (chain B), 274GLN and 350ILE (chain C), 157GLY (chain D), 
186V AL and 459MET (chain E) and 321ALA (chain F).For other chains (B, C, D, E and F) of TA2Fi protein, it is seen in 
the bottom figures that the deviated maxima are: 61GLY, 328ALA and 448ASP (chain B), 57ASN, 214GLN and 466ALA 
(chain C), 105ALA and 267ALA (chain D), 261LEU, 330GLU and 440ASN (chain E), 72ASP, 261LEU and 353ALA 
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(chain F). Some deviated minima are: 457GLU, 474HIS and 492GLU (chain B), 163GLN and 233PRO (chain C), 24PRO, 
409VAL and 450GLY (chain D), 208PHY (chain E), 450GLY (chain F). 
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Figure 2.5: Comparison of theoretical and experimental 
B-factors for the chain B of a subunit of MF! (upper 
figure) and TA2F( (lower figure) proteins. B-factors 
are plotted as a function of residue index. The theoretical 
data were obtained from calculation using oGNM server 
[26] and experimental data were reported in the PDB files 
(ID: lbmf for MF! and ID:2QE7for TA2F[). Sigma Plot 
was used to produce the curves. 
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Fig.2. 7. Comparison of theoretical and experimental 
B-factors for the chain C of a subunit of MF! (upper 
figure) and TA2F[ (lower figure) proteins. B-factors 
are plotted as a function of residue index. The 
theoretical data were obtained from calculation using 
oGNM server [26] and experimental data were reported 
in the PDB files (ID: lbmf for MFi and ID:2QE7 for TA2F0. 
Sigma Plot was used to produce the curves. 
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Figure 2.6: Comparison of theoretical and experimenta 
1 B-factors for the chain B of a subunit of MF[ (upper 
figure) and TA2F[(lower figure) proteins. B-factors 
are plotted as a functionof residue index. The theoretical 
data were obtained from calculation using oGNM server 
[26] and experimental data were reported in the PDB files 
(ID: lbmf for MFi and ID:2QE7 for TA2Fi). Sigma Plot 
was used to produce the curves. 
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Fig.2. 8. Comparison of theoretical and experimental 
B-factors for the chain D of [3 subunit of MFj 
(upper figure) and TA2F[ (lower figure) proteins 
B-factors are plotted as a function of residue index. 
The theoretical data were obtained from calculation 
using oGNM server [26] and experimental data 
were reported in the PDB files (ID: lbmf for MF[ 
and ID:2QE7 for TA2F0. Sigma Plot was used to 
produce the curves. 
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Fig.2. 9. Comparison of theoretical and experimental B-factors 
for the chain E of (3 subunit of MF { (upper figure) and 
TA2Fi (lower figure) proteins. B-factors are plotted as a 
function of residue index. The theoretical data were 
obtained from calculation using oGNM server [26] and 
experimental data were reported in the PDB files (ID:lbmf 
for MF, and ID:2QE7 for TA2F,). Sigma Plot was used to 
produce the curves. 
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Fig.2. lO.Comparison of theoretical and experimental 
B-factors for the chain E of (3 subunit of MF! (upper 
figure) and TA2Fj (lower figure) proteins. B-factors 
are plotted as a function of residue index. The 
theoretical data were obtained from calculation using 
oGNM server [26] and experimental data were 
reported in the PDB files (ID:lbmf for MFj and 
ID:2QE7 forTA2F[). Sigma Plot was used to 
produce the curves. 
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Fig.2. 11. Comparison of theoretical and experimental 
B-factors for the chain D of [3 subunit of MF! (upper 
figure) and TA2F[ (lower figure) proteins. B-factors 
are plotted as a function of residue index. The theoretical 
data were obtained from calculation using oGNM 
server [26] and experimental data were reported in 
the PDB files (ID:lbmf for MF, and ID:2QE7 for 
TA2Fi). Sigma Plot was used to produce the curves. 



Fig.2. 12. Comparison of theoretical and experimental 
B-factors for the chain D of (3 subunit of MF[ (upper 
figure) and TA2Fj (lower figure) proteins. B-factors 
are plotted as a function of residue index. The 
theoretical data were obtained from calculation using 
oGNM server [26] and experimental data were 
reported in the PDB files (ID:lbmf for MFj and 
ID:2QE7 for TA2F,). Sigma Plot was used to 
produce the curves. 



The pattern of B factor curves in Fig.2. 1 1 obtained for the G chain (y subunit) of both the MF , and TA2F[ proteins are 
far different from other curves in Figs. 2. 5-2. 10. In the pink curve drawn from experimental data reveals a minimum around 
8ARG-21LYS, a maximum around 35GLU-44TYR and a minimum around 77LEU-81ILE. The blue curve obtained by 
plotting the data obtained by GNM calculation is almost similar to that of the experimental curve but there are very high 
deviations in maxima, minima or flat regions. Although the experimental and theoretical curve pattern for TA2Fi protein is 
similar but the calculated B factors are much higher than those of experimental ones. Remarkably deviated maxima and 
minima in the G (y subunit) chain are: 65HIS, 130GLY, 162ASP and88ALA, 178PRO, 245MET respectively. Similar to the 
chain G, the theoretical B factors for chain H (e subunit) are very high compared to the experimental values. Some highly 
deviated maxima and minima in this chain are: 30GLY, 1 10THR and 106LYS and 64SER respectively. 

III. DISCUSSION 

In the structure of mitochondrial F r ATPase (MF]) determined in the presence of nucleotide and magnesium, three a 
and three P subunits are alternately arranged to form a 013P3 cylinder in which the subunits are folded identically [7]. An 
antiparallel a-helical coiled coil of amino- (N-) and carboxy- (C-) termini of the ysubunit forms the axle of rotation that 
deeply penetrates the 013P3 cylinder of FpATPase. This axle is held at the top orifice and the bottom of the 013P3 cylinder as 
shown in Fig.2. 13 for MF] protein. The bottom part is made of a hydrophobic "sleeve" that acts as a bearing [7]. The 
structure of 013P3 cylinder of F r ATPase from thermophilicBac/W«.s PS3 [41] was found to be similar to that of MF]. The aand 
P subunits of TA2F] [7] forming the (X3P3 cylinder have a similar overall fold as MF[ and TF) structures [7, 41]. It is thus 
clear that the structures of a 3 p 3 subunits of the three F r ATPases as mentioned above are similar. Each a or P subunit is 
composed of an N-terminal six stranded P-barrel, a central domain containing nucleotide binding sites and a C-terminal 
domain [7]. 

The fluctuations exhibited by different chains of a (A, B, C) and P (D, E, F) subunits of TA2F, protein is very low 
compared to the fluctuations of different chains of the MF[ protein as observed in red and black curves in Figs. 2. 1-2.2. The 
fluctuations hidden in the red curves are revealed in the expanded green curves. In the fluctuation curves of B and E chains 
of TA2Fi protein, very high peaks appeared near 400 residue position. Considering these points, it seems to be that the 
structures of a and P subunits of TA2F t protein are very rigid compared to the MF! protein. The differences in fluctuations of 
these two proteins may be due to their structural differences. The structure of TA2Fi protein was determined in the absence 
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of nucleotides and magnesium so that the a and P subunits could not bind any nucleotide and they remained in open 
conformations. The structure of MFi protein was determined in the presence of nucleotides and magnesium and 
consequently the a and P subunits remained in closed conformation by binding Mg-nucleotide except one P subunit which 
remained empty and open. 




Fig. 2. 13. The chain A of a subunit and the chain F of (3 subunit together with y subunit (chain G) of MFi [26] are shown. 
The chains are represented by letters in the circles. The carboxy (C) terminal part of the y subunit deeply penetrates into 
the molecular bearing made of the bottom part of the a 3 (3 3 subunits. The upper and lower circles represent the entrance 
(orifice) and bearing regions of the a 3 (3 3 cylinder. The figure is made using Jmol software [34]. 

In the case of MFi, the amino terminal regions through the A-F chains of aP subunits around 24-95 residues is almost 
rigid (Figs. 2. 1-2.2). The regions 100-400 residues in each fluctuation curve are more flexible compared to the amino 
terminal region. The regions around 400-500 residues containing carboxyl terminal regions show the highest flexibility 
compared to the other two regions. The curve for the chain B exhibits the highest mobility and the region around 400-500 
residues in the curve for the chain F exhibits less flexibility. The a subunits are more flexible compared to the P subunits. 

In some mutational studies [31, 32], the residues 155LYS and 182ARG in the P subunit of the Escherichia coli¥\- 
ATPase (EFi) corresponding to the residues 162LYS and 189ARG in the P subunit of MF] were found important for 
generation and stabilization of the catalytic transition state. The residue P-182ARG plays role in the required conformational 
change at the catalytic sites for multisite catalysis [34, 35]. The residues 162LYS and 189ARG in the P subunit of MFj are 
located within the minima in the black curves obtained for chains D, E and F in Fig. 2. 2. These minima have some 
fluctuations i.e. flexibility. The residues in the green curve for D, E and F chains of TA2F] protein corresponding to the 
residues 162LYS and 189ARG in MFj also lie in minima but with some fluctuations/flexibility region. These similarities 
indicate that some flexibility/conformational changes are necessary for catalysis in FpATPases irrespective of their origins. 
The residue 188GLU in the P subunit of MFj [19] lies in a position to activate the water molecule that enhances a 
nucleophilic attack on the terminal phosphate. In the case of FpATPase from Bacillus PS3 (TFi), mutations of PE190D and 
PE190DCax residues (corresponding to P188GLU in MFi) severely impaired catalysis [40]. The position of the residue 
P188GLU in the black curves for D, E and F chains of MFi is within minima with some fluctuations. The position of the 
residue in the green curves for TA2Fj protein corresponding to P188GLU in MFi remains in the regions having some 
fluctuations. Thus it is evident that some sorts of molecular motion or change in conformation occur in the relevant sites of 
subunits of F r ATPase for nucleotide hydrolysis. 

Considering the overall fluctuations in the Figs. 2. 1-2.2, it is clear that a and P subunits of a 3 p 3 cylinder in TA2F, 
protein are highly rigid compared to the flexible a 3 p 3 cylinder of MF[ protein. The y (G) subunit of MFj protein almost fully 
penetrates into the cavity of a 3 p 3 cylinder up to the molecular bearing (Fig.2.13) and its motion may be highly restricted by 
the a 3 p 3 cylinder so that it shows very little fluctuations (Fig.2.3) compared to the a 3 p 3 subunits. In the case of TA2Fj 
protein, only a part of the y (G) subunit of TA2Fj protein enter through the orifice and it remains in the upper portion of the 
cavity of the a 3 p 3 cylinder (Fig.2.14). Rest of the cavity remains empty. Moreover, the protruding portion of the y (G) 
subunit is attached to the e (H) subunit. The y (G) subunit of TA2F[ protein thus feels very low level of restrictions from a 3 p 3 
cylinder around it and this may be an important reason for very high fluctuations of the y (G) subunit. The £ (H) subunit 
remains completely outside the cavity of a 3 p 3 cylinder so that it gets no resistance from a 3 p 3 cylinder and it can move more 
freely during molecular motion. This may be one reason why e (H) subunit shows exceptionally high fluctuations compared 
to the all other subunits of TA2F, protein and that of MF] protein. 




Fig.2.14. The chain A of a subunit and the chain F of (3 subunit together with y subunit (chain G) and s subunit (chain H) of 
TA2F! [24] are shown. The chains are represented by letters in the circles. The carboxy (C) terminal part of the y subunit 
remains in the cavity region of the <x 3 [3 3 cylinder. The upper and lower circles represent the entrance (orifice) and bearing 
regions of the <x 3 (S 3 cylinder. The figure is made using Jmol software [34]. 



www.ijmer.com 



3490 I Page 



International Journal of Modern Engineering Research (IJMER) 
www.iimer.com Vol. 3, Issue. 6, Nov - Dec. 2013 pp-3483-3492 ISSN: 2249-6645 



The more localized atoms generate relatively low B factor values and they usually take part to form secondary 
structures, neighboring disulfide bridges, etc. The chain termini exhibit fluctuation peaks in several proteins due to their 
inherent enhanced flexibility. The surface atoms generally exhibit larger fluctuations, compared to core atoms. In the case of 
a multimer protein consisting of several subunits, theoretical curves fall above the experimental curves at the regions of 
subunit-subunit interfaces, due to the absence of interactions between neighboring subunits [9]. 

The y subunit of F r ATPase rotates inside the cylinder made of a 3 p 3 subunits and rotation was first observed in TFj 
[26] using various probes such as actin filament [26], a spherical bead or a bead duplex [30, 31] and so on. Subsequently, 
rotation was also observed in F r ATPases from Escherichia coli (EF])and from chloroplast {CF { ) [32, 33]. Rotational 
experiment on TA2Fj protein has not yet been reported. According to a push pull mechanism [51], when a P subunit binds a 
nucleotide its upper portion bends toward and apparently pushes the axle at the top and when the P subunit releases the 
nucleotide it pulls the axle. The push pull actions generate rotation of the slightly curved axle in a conical fashion. This 
mechanism requires that the axle to be rigid and its tip held in the bearing acts as a pivot. 

In MFj, the sleeve/bearing consists of the six proline-rich loops (a287-294 and P274-281) of the 013P3 cylinder. The 
residues 287-294 remain in the minima regions of black curves obtained for chains A, B and C of a subunit as seen in 
Fig.2.1. Comparing the green curves of chains A, B and C of a subunit of TA2Fi protein in the same figure it is observed 
that the corresponding residues for TA2Fj protein also lie in the minima of green curves. Similarly, the residues 274-281 in 
the P subunit of MF, remain in the minima of black curves and the corresponding residues in TA2F] protein lie almost in the 
minima of green curves as shown in Fig. 2. 2 plotted for the chains D, E and F. Thus from the study of molecular motion it 
appears that the residues forming the molecular bearing of both MF] and TA2F] proteins remain in the minima which implies 
that bearing is quite rigid that helps to hold the carboxy (C) terminal tip of the axle giving an overall structural stability of F r 
ATPases. However, an axle-less mutant of F r ATPase from thermophilic bacterium (TFj) could rotate for >100 revolutions 
[29] which implies that neither a rigid axle nor a fixed support at the bottom is necessary for rotation. 

The highly conserved and charged DELSEED sequence (P394ASP-400ASP) of P subunit that interacts with the 
ysubunit can be replaced with alanine without affecting rotational torque [28] giving an impression that the interactions 
between the protruding portion of y subunit and the top surface of P subunit is less important. As seen in Fig. 2. 2, the residues 
394ASP-400ASP of the P subunit of MFj represented by black lines for the chains D, E and F lie within minima. The 
corresponding residues for the TA2Fj protein also remain in the minima regions of green curves. This suggests that these 
residues are restricted or least mobile. 

IV. CONCLUSION 

Fl-ATPase is a motor protein that generates rotation during ATP hydrolysis. A motor protein is a machine that 
converts chemical energy into mechanical energy and forms the basis of a specific cellular function. The proteins undergo 
motion and exhibit different ranges of flexibility during their functions. The molecular motions of two proteins of different 
sources namely, Fl-ATPase from bovine mitochondria (MF1) and Fl-ATPase from thermoalkaliphilic bacteria (TA2F1) 
were studied in this thesis using Gaussian Network Model (GNM). The data obtained by calculation using GNM server [26] 
for MF1 and TA2F1 proteins for the 1st slowest mode of vibration were compared and analyzed. A force constant y= 0.054 
kcal/(mole-A 2 ) was obtained for TA2F1 protein and y= 0.08 kcal/(mole-A 2 ) [27] was for MF1 protein. The variation in the 
force constant is =0.026 kcal/(mole-A 2 ) which implies that the inter residue contact potential is stronger in MF1 than in 
TA2F1. The degree of fluctuations of the different chains of a subunit (A, B and C) and /? subunit (D, E and F) of the MF1 
and TA2F1 proteins varies from 0~ 2.5A2 (Figs. 2. 1-2.2). The largest fluctuations ~ 2.5 A 2 in MF1 occurs in the carboxy 
terminal region of the chain B. This suggests that a 3 /? 3 the cylinder of MF1 is flexible whereas the a 3 /? 3 cylinder of TA2F1 is 
highly rigid. The y subunit of TA2F1 remains in the cavity of the a 3 /? 3 cylinder near the orifice (Fig. 2. 14) and it gets less 
resistance during motion compared to the y (G) subunit of MF1 (Fig. 2. 13). This may cause a fluctuation ~ 2.5A 2 higher in 
TA2F1 than in MF1. The fluctuation of £ (H) subunit of TA2F1 tagged with the protruding portion of the y (G) subunit is the 
highest -3.5 A 2 compared to all other subunits of both the proteins. The overall agreement between theoretical and 
experimental values of B factors is good which is supported by the correlation coefficient. The calculated correlation 
coefficient for MF1 was 0.71 (71%) [27] and for TA2F1 was 0.64 (64%). The residue /M88GLU in MF1 responsible for 
ATP hydrolysis andthe corresponding residue in TA2F1 exhibit flexibility. The hydrophobic 'sleeve '44 forming the 
molecular bearing was found to be located in the minima for both the proteins. 
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